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Abstract 

Protein mobility is conventionally analyzed in terms of an effective diffu¬ 
sion. Yet, this description often fails to properly distinguish and evaluate 
the physical parameters (such as the membrane friction) and the biochemi¬ 
cal interactions governing the motion. Here, we present a method combining 
high-density single-molecule imaging and statistical inference to separately 
map the diffusion and energy landscapes of membrane proteins across the 
cell surface at ~ 100 nm resolution (with acquisition of a few minutes). 
When applying these analytical tools to glycine neurotransmitter receptors 
(GlyRs) at inhibitory synapses, we find that gephyrin scaffolds act as shallow 
energy traps (~3 fcgT) for GlyRs, with a depth modulated by the biochemi¬ 
cal properties of the receptor-gephyrin interaction loop. In turn, the inferred 
maps can be used to simulate the dynamics of proteins in the membrane, 
from the level of individual receptors to that of the population, and thereby, 
to model the stochastic fluctuations of physiological parameters (such as the 
number of receptors at synapses). Overall, our approach provides a powerful 
and comprehensive framework with which to analyze biochemical interac¬ 
tions in living cells and to decipher the multi-scale dynamics of biomolecules 
in complex cellular environments. 

Key words: Single Molecule, Neurons, Statistical Physics 
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Introduction 

Determining the parameters that regulate the mobility of proteins in cells is 
key for many cellular functions. The motion of proteins depends on a variety 
of factors, including the local viscosity, their intermittent binding to other 
proteins, the molecular crowding and the dimensionality of the accessible 
space (pj) . Since all these factors are difficult or impossible to reconstitute 
in vitro using purified constituents, there is a compelling need for analytical 
tools that bypass in vitro assays and directly access the properties of macro- 
molecular assemblies and the kinetics of their interactions in their native 
cellular environment. 

Thanks to single-molecule (SM) imaging tools, it is now possible to 
record trajectories of individual proteins in a variety of cellular systems. 
An important challenge is to extract relevant biochemical and biophysical 
information from these trajectories. This is commonly done by computing 
the mean square displacement (MSD) along the trajectories and estimating 
the effective diffusion coefficient of the molecule. By associating the diffu- 
sional states to the functional states of the biomolecules, one can identify 
molecular behaviors (1) and evaluate the transition kinetics between them 
(2). Although this approach has often proved useful, it is conceptually in¬ 
appropriate in many biological situations. Measuring a diffusion coefficient 
places emphasis on the friction encountered by the protein and assumes that 
the movement is characterized by a MSD scaling linearly with time. Yet, 
the primary factor controlling the motion of a protein is often not the fric¬ 
tion but, rather, its interactions with molecular or macromolecular partners 
leading to transient stabilization or transport. In this case, the relevant in¬ 
formation is not the diffusion coefficient but the binding energies between 
the protein of interest and its interacting partners. Furthermore, regulatory 
processes are often mediated by changes in these binding energies, which 
should ideally be evaluated with in situ measurements. 

Methods that go beyond the computation of the MSD generally aim to 
identify deviations from Brownian movement within single molecule trajec¬ 
tories, due for instance to trapping or transport (0-0). However these meth¬ 
ods essentially remain ad hoc tools and do not constitute a comprehensive 
framework to describe the parameters underlying the motion. Furthermore, 
biological media are often spatially inhomogeneous and this heterogeneity 
is poorly conveyed by measuring a few, sparse trajectories. A conceptu¬ 
ally different approach using Bayesian inference methods has been recently 
proposed to analyze the motion of molecules 6,0). It assumes that the 
membrane environment is characterized by two spatially-varying quantities: 
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(i) the diffusivity D (r) = ksT/^igc) (where 7 (r) is the local viscosity), (ii) 
the potential energy V (r) that reflects the biochemical interactions of the 
molecule. In this framework, the protein is a random walker with a motion 
governed by the Langevin equation 0 ): 


dr 

dt 


D (r) VI/ (r) 

k^f 


+ y 7 2D (r)£ (t) 


(1) 


(where £ (t) is a rapidly varying Gaussian noise with zero mean). From 
a general standpoint, a knowledge of D (r) and V (r), which are protein- 
specific, can not only reveal how fast the protein moves in the membrane 
but also identify areas where it can be stabilized (energy traps) or from 
which it is excluded (energy barriers). However, in the few cases where 
D (r) and V (r) have been experimentally determined 4M), the analysis 
has been limited to movements confined in local regions (< 1 / 1 m 2 ), falling 
short of providing a complete description of the heterogeneous diffusivity 
and energy landscapes in the cell membrane. 

Here, we introduce a novel and generic approach, combining high density 
single molecule imaging and computational tools, that enables the mapping 
of the environment of membrane receptors across the entire cell surface and 
at ~ 100 nrn resolution. This approach allows the mapping of the membrane 
over regions of several hundred pm 2 in a few minutes of data acquisition. 
Furthermore, the inferred maps are used to numerically generate massive 
number of trajectories. These simulated trajectories, whose characteristics 
match those of the experimental ones, enable a complete analysis of the 
dynamics in the complex membrane environment by means of various sta¬ 
tistical estimators. To illustrate the relevance and benefits of our approach, 
we applied it to the neuronal membrane, a cellular system in which the 
spatial organization is critical for the detection and processing of external 
information. In past years, tracking experiments have underlined the role 
of membrane dynamics in ensuring rapid exchange of receptors (e.g. glu¬ 
tamate, glycine or GABA receptors) between extrasynaptic and synaptic 
localizations 0 )- Therefore, the number of receptors at synapses depends 
on the motion of receptors at the cell surface and their stabilization at synap¬ 
tic loci, the latter being regulated by the number of scaffolding molecules 
and the affinity of the receptor-scaffold interactions 0). A quantitative 
analysis of the protein mobilities and of their regulatory mechanisms is thus 
paramount for characterizing and modeling the variability of the synaptic 
response and the plasticity of the nervous system (involved in higher brain 
functions such as learning and memory or during pathological processes). 
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Results and discussion 


Mapping the diffusion and energy landscapes with Bayesian infer¬ 
ence. Our approach for the large-scale mapping of D (r) and V (r) builds 
on Bayesian statistical tools recently developed to analyze the motion of 
individual particles 0,0). The principle of the method is as follows (see de¬ 
tails in the Supplementary materials). We first acquire high-density single- 
molecule data (14, HE), with a number of individual translocations on the 
order of 1000 — 10000 /pm 2 . Next, the surface of the cell is meshed with 
sub-domains S t ,j (labeled with the index (i,j) along the x and y axis) with 
a size proportional by a factor 6 ~ 2-3 to the average step size of a translo¬ 
cation, such that consecutive positions of the molecules are either in the 
same or in adjacent domains (Fig. la). From the information contained 
in the massive number of individual translocations, we determine D,j and 
VVjj in each sub-domain (i,j) using Bayesian inference techniques adapted 
from 0 ). In brief, we compute the global posterior distribution P of the pa¬ 
rameters and {VV(;j}(,; given the observed trajectories 

Since all the sub-domains are independent, P is the product of the posterior 
distributions inside each of them: 


P 


({VVi J -} (ij .),{A J -} (ij - ) |{T fc } (fc) )= ) xP(VV,D) (2) 
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(3) 




where p designates the index for which the points rjj of the k th trajectory 
are in S t j, o is the experimental localization accuracy (~30 nm), At the 
acquisition time and P (VV, D) the prior information on the potential and 
the diffusivities. In the second line of equation (2) we display the prior 
we commonly used, Jefferey’s prior, that is discussed in the supplementary 
Materials. The estimators (Z?^ AP , VV^ I 4 AP ) of the local diffusivity and 
force are the Maximum a Posteriori (MAP) of the posterior distribution P 
0,0. Finally, we solve the inverse problem to determine in each sub- 
domain the potential field Vij associated to the force. The estimation of Vi j 
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Figure 1: General scheme of the assay, a) Principle of the Bayesian in¬ 
ference method. Left: High-density single-molecule data (red dots) are 
recorded at the cell surface. Right: In a mesh domain, multiple translo¬ 
cations (top) are used to infer the local diffusivity and force (gradient of 
the potential) that underlie the motion (bottom), b) GlyRs (blue) diffuse 
in the membrane and are in dynamic equilibrium between synaptic and ex- 
trasynaptic domains in the neuronal membrane. At synapses, GlyRs are 
stabilized by their interactions with gephyrin clusters (orange), which can 
be modeled as trapping potential (with depth Us)- c) Expression constructs 
of transmembrane (TM) proteins with an extracellular pHluorin tag and an 
intracellular interaction loop derived from the GlyR /3-subunit, d) Principle 
of high-density single-molecule uPaint imaging (|15l ). 
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is performed by minimizing £({Vi,j}), defined as: 

£ m,j}\ (hi) G {N {hi)} ^ °) = E ( W m - vv*f AP ) 2 +p ( 5 ) £ (W^) 2 

(*,i) (*,j) 

(4) 

with iV (/, j) the number of neighboring occupied mesh domains, (3 (<5) a con¬ 
stant (optimized on numerically generated trajectories) depending on <5 (see 
Supplementary Information). Eventually, the set of quantities {E>^ AP , VV)j} 
constitute the diffusivity and potential energy maps. 

Glycine receptors and their interactions with scaffolding proteins. 

We applied our inference-based mapping method to investigate the dynam¬ 
ics of GlyRs in the neuronal membrane as well as their stabilization at 
inhibitory synapses 0 - This stabilization is achieved through the binding 
of the receptors to the scaffold protein gephyrin (Fig. lb) via an intracellular 
loop (the /3-loop) present in the two /3-subunits of the pentameric GlyR com¬ 
plex. The high affinity component of the /3-loop-gephyrin interaction is in 
the nanomolar range {Kjj ~ 20 nM), as determined by isothermal titration 
calorimetry 0). To characterize the GlyR-gephyrin interaction in living 
neurons, we used recombinant membrane proteins consisting of a trans¬ 
membrane domain (TM), a C-terminal pHluorin tag (a pH-sensitive GFP 
mutant that is quenched in intracellular acidic vesicular compartments), 
that were fused N-terminally to the intracellular GlyR /3-loop (Fig. lc). 
This /3 WT -TM-pHluorin construct recapitulates the interactions of the en¬ 
dogenous GlyR complexes with the gephyrin scaffold proteins, with the im¬ 
portant benefit that individual elements of the receptor-scaffold interaction 
can be manipulated independently 0). It also overcomes the difficulty of 
defining the sub-unit composition of oligomeric receptors where transfected 
sub-units compete with endogenous ones. As a control, we used /3 - -TM- 
pHluorin, a construct with a mutated /3-loop that does not interact with 
gephyrin. 


High-density single-molecule imaging of TM proteins. We acquired 
a high-density of individual trajectories using uPaint, a single-molecule tech¬ 
nique in which cells are imaged at an oblique illumination in a buffer con¬ 
taining dye-labeled primary antibodies (1151 ). As antibodies (in our case, 
anti-GFP antibodies coupled to Atto647N dyes) continuously bind to their 
membrane targets, they can be tracked until they either dissociate or photo¬ 
bleach (Fig. Id and Supplementary Videos 1 and 2). Hence, the entire field 
of view is constantly replenished with new fluorescent labels and a large 
number of individual trajectories covering a field of view of ~ 500-1000 
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pm? can be recorded. Experiments were performed on cultured rat hip¬ 
pocampal neurons co-transfected with mRFP-tagged gephyrin and with the 
pHluorin-tagged transmembrane constructs (Fig. lc). In typical measure¬ 
ments, movies were recorded for ~5-15 minutes with an acquisition time At 
= 50 ms (Supplementary Videos 1 and 2), yielding up to hundreds of thou¬ 
sands of individual translocations per field of view, with an average of 30 
points per mesh domain (size ~100xl00 nm 2 ). On this time scale, the cells 
and synaptic sites remained relatively stable, meaning that the diffusivity 
and energy landscapes could be considered constant. 

Diffusion and energy maps of TM proteins. Figures 2a-f show ex¬ 
amples of the diffusivity and energy maps for the two constructs /3 wt -TM- 
pHluorin and /3~-TM-pHluorin. In both cases, the diffusion map exhibits 
fluctuations at short scale (~1 pm or less), with local peaks and valleys and 
a characteristic diffusivity in the range of 0.05-0.2 ^m’.s' 1 (Fig. 2b and 2e). 
More striking differences were observed between the energy landscapes. For 
/3 WT _TM, the landscape is characterized by the existence of small regions 
(< 0.5 /mi 2 ) corresponding to local energy minima (Fig. 2c). Importantly, 
gephyrin clusters coincide with energy minima, consistent with the stabiliza¬ 
tion of the transmembrane proteins at synaptic sites. Yet, we also observed 
that some other minima did not colocalize with gephyrin clusters, suggest¬ 
ing that /3 WT -TM-pHluorin might interact with other partners outside of 
synapses (such as the cytoskeleton or lipid domains). It is possible that 
these extrasynaptic interactions are still mediated by gephyrin (present in 
number too small to be detected), since gephyrin is known to associate with 
GlyRs both inside and outside of synapses 0). In contrast, the energy 
map for /3 _ -TM (Fig. 2f) shows variations at a longer length-scale, without 
correlation to gephyrin clusters. 

To more quantitatively compare the heterogeneous properties of the neu¬ 
ronal membrane for /3 wt -TM and /3 - -TM, we computed two quantities (av¬ 
eraged over 7 cells in each case): (i) the distribution of diffusion coefficients 
in the maps (Fig. 2g), and (ii) the rugosity of the energy landscape, defined 
as the standard deviation of the potential in circular region of defined radius 
averaged over the complete surface of the cell over circular region of the fluc¬ 
tuations of the potential averaged over circular regions of increasing radius 
(Fig. 2h and Supplementary Information). These parameters revealed that 
the interacting /3-loop led to a lower average diffusivity (0.06 /mn 2 .s _1 and 
0.13 /xm 2 .s -1 for /3 wt -TM and /3 - -TM, respectively) and a larger rugosity of 
the potential. This is consistent with the notion that moving TM proteins, 
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Figure 2: Diffusion and energy maps in live neurons, a) Fluorescence im¬ 
ages of cultured neurons expressing mRFP-gephyrin and /3 WT -TM-pHluorin. 
Scale bar: 10 /im. b-c) Diffusion and energy maps, d-f) Equivalent set of 
images and maps for /3”-TM-pHluorin. g) Distribution of diffusion coef¬ 
ficients for the membrane constructs /3 wt -TM (black), ^ S ' 403D -TM (blue) 
and /3 - -TM (red). The vertical bars on the x-axis indicate the mean values 
of the respective distributions. Insert: distribution in a lin-log scale, h) 
Rugosity of the membrane potential as a function of the region radius. 
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when bound to intracellular scaffolding proteins, encounter more obstacles 
that increase the viscosity of their environment. Also they are more likely 
to interact with membrane or sub-membrane structures that contribute to 
the energy landscapes. 


Synaptic scaffolds as crowded energy traps. Given the pronounced 
differences between the energy landscapes of the /3 WI -TM and /3 - -TM con¬ 
structs, we examined the behavior of /3 wt -TM at gephyrin clusters in closer 
details. An example of the energy profile of /3 wt -TM proteins at a synaptic 
cluster (identified by the presence of mRFP-gephyrin fluorescence) is shown 
in Fig 3a. The profile reinforces the view that clusters of scaffolding proteins 
act as energy traps for membrane receptors (fill . 12, 18). The average trap 
depth was 3.6 ± 0.4 ksT (± S.E.M., n = 69 clusters), a relatively shal¬ 
low potential from which receptors can escape rapidly. Yet, about 15% of 
clusters had stabilization energies greater than 6 ksT, corresponding to a 
much more stable anchoring of receptors (Fig. 3b). This reflects the hetero¬ 
geneity of the synaptic domains in the neuronal membrane and underlines 
the need for measurements at the single synapse level. Of note, the binding 
energies between /3 wt -TM and gephyrin seem to be significantly lower than 
the stabilization energy of AMPA receptors at synaptic sites, for which 25 
% of the wells had a depth larger than 8 TTbT (Hi!)- The method used in 
®), also based on a combination of high-density single molecule imaging 
and statistical inference, evaluates the diffusion and drift by computing the 
maximal likelihood estimation in a mesh square as described in ( 3 ). The 
confining potentials were subsequently evaluated by L2 minimization of a 
parabolic shaped potential from the force (drift) fields. In (Ell ) the authors 
do not discuss the role of known biases with confining potentials (see JM)) 
or the effect of the positioning noise, and do not provide information on the 
posterior distribution of the parameters. It is thus delicate to precisely com¬ 
pare their experimental results with ours. Yet, given that the diffusivity of 
AMPARs at excitatory synapses appears to be higher than the diffusivity of 
GlyRs at inhibitory synapses (gephyrin clusters), higher confining potentials 
may be necessary to stabilize the AMPARs. In addition, we noticed that 
the average diffusivity of /3 wt -TM (~0.01 pm 2 s~ 1 ) inside gephyrin clusters 
was reduced by a factor ~6 compared to extrasynaptic regions (Fig. 3c), 
probably due to the combined effect of membrane crowding within synaptic 
sites and the binding to scaffolding elements. In comparison, the diffusivity 
of /3“-TM proteins inside gephyrin clusters, which we expect to be predom¬ 
inantly influenced by molecular crowding (1211 ). was 0.07 /im 2 .s _1 (Fig. 3c), 
only a factor ~2 lower than in extrasynaptic domains. In other words, the 
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Figure 3: Analysis of the synaptic gephyrin scaffold, a) Example of a 
gephyrin cluster (indicated by a box) acting as a local trap in the energy 
landscape. Scale bar: 5 pm b) Cumulative distribution function (c.d.f) of 
trapping energy for the constructs /3 wt -TM (black) and /3 S403£) -TM (blue). 
Insert: mean values of the distribution. Error bars indicate the S.E.M. c) 
Mean diffusivity for /3 wt -TM (black), /3' S403J:) -TM (blue) and /3 _ -TM (red). 
Error bars indicate the S.E.M. 


synaptic scaffold stabilizes the receptor by simultaneously diminishing the 
diffusivity of the receptor and by acting as a trapping potential. 

Modulation of the /3-loop gephyrin binding affinity. Since the compu¬ 
tation of the energy landscape allows the unambiguous distinction between 
interacting membrane constructs and those lacking interaction domains, we 
tested the sensitivity of our approach with the phosphomimetic construct 
^S'403 D_jm, a mutated /3-loop known to have a lower gephyrin binding affin¬ 
ity in vitro (Kd ~ 0.9 pM (19)) (Fig. lc). As a result, /3‘ S ' 403D -TM displayed 
increased membrane diffusion and reduced synaptic accumulation compared 
to /3 wt -TM. The phosphorylation of amino acid residue S403 of the GlyR/3 
subunit by protein kinase C thus contributes to the regulation of GlyR levels 
at inhibitory synapses 0). The diffusion and energy landscapes of {3 S403D ~ 
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TM (computed over 6 different cells) yielded a diffusivity (average value 0.11 
/rm 2 .s _1 ) and an energetic rugosity precisely intermediate between those of 
the wild-type and of the binding-deficient constructs (Fig. 2g-h). Com¬ 
pared to /3 wt -TM, the average trap depth of ^‘ S403£) -TM at synaptic sites 
was reduced to 2.4 ± 0.4 ksT (n = 58 clusters), with less than 5% of the 
traps above 6 ksT (Fig. 3b). Inside clusters, the average diffusivity (0.015 
^m 2 .s _1 ) was slightly higher than for the wild-type (Fig. 3c). 

Importantly, the binding energy reported here corresponds to TM pro¬ 
teins moving in a two-dimensional membrane and interacting with macro- 
molecular gephyrin scaffolds that are believed to be two-dimensional as well 
( 23 , [ 23 I ). This is in contrast with measurement of the equilibrium con¬ 
stant Kjo by isothermal calorimetry, which reports on the individual inter¬ 
action between the /3-loop and the scaffolding protein in an isotropic, three- 
dimensional measurement of the /3-loop-scaffold interaction. Obtaining the 
stabilization energy thus constitutes a first and important step to bridge the 
gap between in vitro and in situ biochemical measurements. When further 
complemented with data on the ultra-structure and stoichiometry of synap¬ 
tic scaffolds (that are now accessible with single molecule imaging techniques 
(0, [jib), we expect our approach to enable a true determination of the two- 
dimensional affinity of the membrane proteins for the synaptic scaffolds 


Connecting the landscapes and the global mobility of proteins. An 

important question for the dynamics of proteins is how the variability of 
their diffusion and energy landscapes at short scale (~ 100 nm) affects their 
long-distance mobility and, thereby, the kinetics of many intermolecular re¬ 
actions. Reaching a multi-scale description of the motion in the membrane 
has long been a challenge in single-molecule experiments. High-density sam¬ 
pling is usually achieved with poorly stable probes, yielding numerous but 
short trajectories 0,0. In contrast, long trajectories obtained with more 
stable markers such as quantum dots ( 261 ) only provide a sparse sampling 
of the cell surface. Furthermore, the nature of the motion, such as sub¬ 
diffusion, may prevent efficient space sampling with single long trajectories. 
Here, we adopted a different strategy and used the inferred maps as phe¬ 
nomenological templates to simulate the motion of proteins. Practically, we 
used the Gillespie scheme (1271 1 to generate individual trajectories lasting up 
to 500s (see Methods and Supplementary Information). 

From a large number of simulated trajectories, we could compute ensemble- 
averaged quantities. We first evaluated the propagator namely the 

probability density function of moving a distance d in a time t, which is 
the fundamental estimator characterizing the random motion in a complex 
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Figure 4: Analysis of simulated trajectories in the inferred maps. Unless 
otherwise mentioned, the results correspond to the constructs /3 wt -TM 
(black), /3 s ' 403£> -TM (blue) and /3~-TM (red), a) Ensemble-averaged propa¬ 
gator II (d,t), defined as the probability density function (pdf) to move by a 
given distance in f=10 s. b) Propagator II(d, t) for the construct /3 wt -TM, 
computed at different times t. The plain lines are adjustments with the gaus- 
sian curves exp(— d 2 /2y 2 (t))/2 itx 2 {t)■ c) Curves x(f). Inset: Propagators 
for the construct /3 wt -TM as a function of the rescaled variable p = d/x(t)- 
d) Mean-squared displacement (MSD) as a function of time. The straight 
lines indicate the sub-diffusive behavior at short time-scales, e) Time-course 
of the number of receptors at a single synapse. Inset: distribution of the 
minimum (in red) and maximum (in blue) number of receptors computed 
over traces of 300 s for all the gephyrin clusters, f) Autocorrelation func¬ 
tions (in grey) for the time traces of number of receptors at gephyrin clusters 
(computed over 300 s). The red line indicates the average autocorrelation 
function. 
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environment (j28j). Although the difference in the average trapping energy 
at gephyrin clusters was only ~1 ksT between /3 wt -TM and /3 S403D -TM, it 
led to significant changes in the mobility, reducing the probability of moving 
over long distances with increasing strength of the /3-loop-gephyrin interac¬ 
tion (Fig. 4a). To more carefully examine the nature of the movement of 
/jWT.'py/^ we plotted n(d, t ) at different times t. The curves could be approx¬ 
imated by gaussian curves exp(—d 2 /2y 2 (t)) with y(t) oc t a and a = 0.33, 
less than 0.5 the value expected for a standard Brownian motion (Fig. 4b-c). 
In fact, these results are consistent with a subdiffusive movement resulting 
from a fractional brownian motion due to heteregoneities in the diffusion 
and energy landscapes (0). Similar results were obtained for ( d s ' 403£ ’-TM 
and /3 - -TM, with a increasing to 0.39 and 0.41, respectively (Fig. 4c). The 
subdiffusive nature of the motion could be further illustrated by computing 
the ensemble-averaged MSD for the three transmembrane constructs (Fig. 
4d). On the time scales 0.05-5 s, all the MSDs increased sublinearly, with an 
anomalous exponent a equal to 0.75, 0.82 and 0.89 for /3 wt -TM, ( g S403D -TM 
and /3”-TM, respectively. The MSD anomalous exponents are slightly larger 
than 2a, likely due to boundary effects associated to the size and geometry 
of the neurons. 

Finally, we examined the implications of the local properties of the mobil¬ 
ity of individual GlyRs on their global distribution in the membrane and on 
the receptor occupancy at synapses. To do so, we simulated the membrane 
dynamics of a population of receptors, using surface densities derived from 
prior experimental reports (Supplementary information). We computed in 
particular the time course of the number of receptors at individual synaptic 
clusters, which we expect to fluctuate due to the exit and entry fluxes of 
receptors (Fig. 4e and Supplementary Information). The exit kinetics at 
a given synapse is determined by the shape and amplitude of the trapping 
potential combined with the reduced diffusivity in the cluster. In contrast, 
the entry kinetics depends on the motion of all the receptors over the entire 
cell surface and need to be computed using the diffusion and energy maps. 
The number of receptors varied significantly over times, as illustrated by the 
distribution of their minimal and maximal numbers at individual synapses 
(Fig. 4e). Furthermore, the time scale of these fluctuations, analyzed by 
computing the autocorrelation function, is comprised between ~ 1 s and a 
few tens of seconds, showing a large heterogeneity among gephyrin clusters 
(Fig. 4f). These observations may account for the dynamic range of receptor 
numbers at synapses and for the variability of synaptic transmission (1 3 il l. 
The receptor fluctuations, that are equivalent to a noise, may also favor the 
transition from one steady state to another during synaptic plasticity 0). 
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Conclusion 

The motion of proteins in the plasma membrane is influenced by both a 
viscous landscape, 7 (r), and an interaction potential, V (r). We have in¬ 
troduced a method to map the interaction energy and diffusion landscapes 
in the cellular membrane with ~ 100 nm resolution over surfaces of several 
hundred pm 2 . The possibility of simulating trajectories in the inferred maps 
offers many possibilities to address the multiscale dynamics of membrane 
proteins. In particular, it bridges the gap between the information obtained 
from numerous, dense but short - trajectories acquired using uPaint or 
sptPALM techniques, and that from the much longer, but usually sparse, 
trajectories extracted through the tracking of proteins labeled with photo¬ 
stable fluorophores (Qdots, nanoparticles). These trajectories can be used to 
accurately evaluate various statistical estimators, thus enabling the analysis 
of the dynamics of biomolecules in complex media. We anticipate that our 
method will be instrumental to identify the factors governing the mobility 
of specific molecules (such as friction, molecular interactions and geometry 
of the cell) and thereby to model and analyze reaction-diffusion processes 
in biological media. As illustrated in the case of GlyR-gephyrin binding, it 
also paves the way to in situ biochemical measurements, which is key for a 
quantitative analysis of the regulation of molecular interactions in a cellular 
environment. Our approach should also be helpful to describe the molec¬ 
ular noise that results from variability of protein concentrations across the 

1 ry an important role in information processing at 



Beyond the case of receptor-scaffold interactions, 


our analytical tools can be applied to other biological questions, such as the 
stability of macromolecular assemblies in the cytoplasm or the nucleus, or 
to the sequence-dependent movement of proteins along DNA (33). 

Materials and methods 

Antibody coupling. Rat anti-GFP monoclonal antibody (Roche) was la¬ 
beled with Atto-647 dye using standard conjugation methods. In brief, 40 
pL of antibodies at 0.4 mg/ml in PBS were mixed with 4 pL of 1 M sodium 
bicarbonate buffer at pH 8.5. This solution was incubated with 10-fold 
molar excess of Atto-647-NHS-ester (Sigma) diluted at 1 mg/rnL in anhy¬ 
drous DMSO. After 1 h of incubation at room temperature, the solution was 
filtered with a Microspin G50 column (GE Healthcare) to remove unconju¬ 
gated dye. The overall coupling efficiency of the dye, estimated by UV-Vis 
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absorption, was about 12%. The labeled antibodies were washed with PBS 
and concentrated using three rounds of centrifugation with a vivaspin500 
10 kDa cut-off PES membrane filter (GE Healthcare). The concentrated 
antibody solution was stored at 4 C and used for up to one week. 

Cell culture and plasmid transfection. Hippocampal neurons from 
Sprague Dawley rats at embryonic day 18 were cultured at a density of 6 
x 10 4 cells/cm 2 on 18 mm coverslips precoated with 80 mg/ml poly-D,L- 
ornithine (Sigma) and 5% fetal calf serum (Invitrogen) as described pre¬ 
viously (1191 1. Cultures were maintained in serum-free Neurobasal medium 
supplemented with lx B27 and 2 mM glutamine (Invitrogen). Cells were 
transfected after 6 to 8 days in vitro using Lipofectamine 2000 (Invitrogen), 
and imaged 1 to 2 days after transfection. All coverslips were co-transfected 
with mRFP-tagged gephyrin and pHluorin-tagged TM constructs, using 0.4 
pg of each plasmid per coverslip. The expression constructs /3 wt -TM- 
pHluorin, ^' S403£l -TM-pHluorin and /3 _ -TM-pHluorin are all described in 
(1191 ). In brief, fj SA0:iD corresponds to the mutation of serine S403 of the 
GlyR/3 subunit that mimics the phosphorylation of the residue by protein 
kinase C. /3 _ -TM corresponds to the double mutation F398A and I400A of 
the wild-type GlyR /3-loop that abolishes binding to gephyrin. 

Cell labeling. Prior to imaging, we prepared a stock solution of diluted an¬ 
tibodies using casein (Vector laboratories) as a blocking reagent. We added 
2 pi of Atto-647 conjugated anti-GFP antibodies and 10 /d of 10 mg/ml 
casein to 40 pi of PBS, resulting in an antibody solution of 0.1-0.2 pM. We 
also prepared a stock of Tetraspeck fluorescent microbeads (Invitrogen) by 
mixing 1 pi of 0.1 pM. microbeads with 400 pi of imaging solution. These 
multi-color fluorescent beads were used as a reference to align the different 
imaging channels and to correct for x/y drifts of the stage and the coverslip. 
The coverslip was mounted in an imaging chamber and incubated with 20 
pi of warmed microbead solution for 10 seconds. After rinsing, the chamber 
was filled with 600 pi of warmed imaging solution (MEMair: phenol red- 
free MEM, glucose 33 mM, HEPES 20 mM, glutamine 2 mM, Na-pyruvate 
1 mM, and B27 lx) and placed on the microscope. To avoid saturating the 
cell membrane with fluorescent antibodies, we first selected a transfected 
neuron and added the fluorescent antibodies at a final concentration of 0.3- 
0.6 nM directly before the start of the acquisition. 

Imaging. Measurements were performed on an inverted epi-fluorescence 
microscope (Olympus 1X70) equipped with a lOOx 1.45NA oil objective and 
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a back-illuminated electron-multiplying CCD camera (Quantem, Roper). 
We imaged the neurons at 37 C in MEMair recording medium using a 
heated stage. For each neuron, we first recorded images of the pHluorin 
signal of the TM constructs and of mRFP-gephyrin fluorescence, using a 
UV lamp (Uvico, Rapp Optoelectronic) and standard sets of filters for GFP 
(excitation 475AF40, dichroic 515DRLP and emission 535AF45) and RFP 
(excitation 580DF30, dichroic 600DRLP and emission 620DF30). Next we 
acquired a uPaint movie of the transmembrane proteins labeled with Atto- 
647-coupled anti-GFP antibodies (20 000 images at 20 frames/s). Atto-647 
dyes were excited with a 640 nm laser and their fluorescence was collected 
through using a 650DRLP dichroic and a 690DF40 emission filter. The laser 
was tightly focused on the back focal plane of the objective. The angle of 
incidence of the beam on the coverslip, controlled by laterally moving the 
focused spot, was just under the limit of total internal reflection, such that 
the laser beam in the sample was almost parallel to the glass surface. This 
angle was slightly adjusted in each experiment to maximize the signal/noise 
ratio of the single fluorescent spots diffusing in the membrane. 


Data analysis. Tracking analysis of the movies was carried out using an 
adapted version of the multiple target tracking algorithm (1341 ). In brief, 
fluorescence spots corresponding to the point-spread function of single emit¬ 
ting fluorophores were fitted with a two-dimensional Gaussian. The centre 
of the fit yielded the position of single molecules with localization accuracy 
~ 30 nm. Trajectories were then computed from individual detections with 
a nearest-neighbor algorithm. 


Simulations in the landscapes. The maps of the diffusion and energy 
landscapes, D (r) and V (r), can be used to simulate the behavior of the 
molecules at different time and space scales. In each mesh sub-domain (i, j ) a 
diffusivity Dij is associated with a potential energy value Vjj. The dynamics 
of the molecules are described by the Fokker-Planck equation: 


dP (r, t|r 0 , t) 
dt 



7 ( r ) 




(5) 

where P(r, t|ro,t) is the conditional transition probability from (ro,to) to 
(r, t). Fokker-Planck equations can always be approximated by Master equa¬ 
tions: 


dP(i,j) ( t ) 
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with in our case 
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(7) 


if the transition is in the x direction and a similar formula in the y direction 
and with W(i.j),(i',j') the transition rate from the ( i',j') site to the (i,j), 
Ax ((Ay)) the mesh size in the x (y) direction, and ^ the poten¬ 
tial gradient acting on the random walker in the x direction when moving 
from (i',j') to (i,j). The motion of the molecule following equation [6] was 
simulated using the Gillespie scheme (j27l ). When the molecule was at the 
site (i,j), the transitions rates, rewritten a u to match Gillespie formalism,;/ 
taking values from 1 to 4, were evaluated on all neighboring sites. We define 
ao = The time, r, to move from the site (i,j) to a neighboring site 

is extracted from an exponential probability density function of rate do, so 
1 ' 1 ' with r\ a random number in [0,1]. The destination 


that r = — log 

ao ° 


site, k, is chosen to satisfy : J2t=o a v — r 2 fl o £ 0 av w ith r 2 a random 

number in [0,1]. Limits of the neuronal cells and unvisited sites are defined 
as inaccessible sites. Note that the trajectory generation process leads to 
trajectories with non-constant time steps. In order to evaluate the different 
estimators, trajectories were regularized to obtain the molecule position at 
regular time lags by imposing that as long as each r was not reached the 
molecule did not move. 
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